Some implications of a new definition of 
the exponential function on time scales 

Jan L. Cieslinski* 

Uniwcrsytct w Bialymstoku, Wydzial Fizyki, ul. Lipowa 41, 15-424 Bialystok, Poland 



Abstract 

We present a new approach to exponential functions on time scales and to 
timescale analogues of ordinary differential equations. We describe in detail 
the Cayley-exponential function and associated trigonometric and hyperbolic 
functions. We show that the Cayley-exponential is related to implicit midpoint 
and trapezoidal rules, similarly as delta and nabla exponential functions are 
related to Euler numerical schemes. Extending these results on any Pade 
approximants, we obtain Pade-exponential functions. Moreover, the exact 
exponential function on time scales is defined. Finally, we present applications 
of the Cayley-exponential function in the g-calculus and suggest a general 
approach to dynamic systems on Lie groups. 
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1 Introduction 

Time scales were introduced in order to unify differential and difference calculus 
[17, 18]. A time scale T is any non-empty closed subset of K. In this paper we con- 
sider the problem of finding good analogues of exponential, hyperbolic and trigono- 
metric functions on time scales, and also finding T-analogues of ordinary differential 
equations (ODEs). In the first part of this paper we present results obtained in [10]. 
Then, many new developments are briefly described. In particular, we show how 
to generate T-analogues of ODEs using known numerical schemes, including the 
implicit midpoint rule, the trapezoidal rule, the discrete gradient method and Pade 
approximants. We introduce exact T-analogues of elementary and special functions 
and exact T-analogues of ODEs. We also present some new results in the g-calculus, 
see [11]. Finally, we show how to define T-analogues of functions on Lie groups using 
the Cayley transformation. 

* e-mail: janekOalpha.uwb.edu.pl 



1 



The classical Cayley transformation, z — > cay(z,a) = (1 — az)~ x (\ + az), is 
a conformal transformation of the complex plane. Generalizations of the Cayley 
transformation on Lie groups and operator spaces are well known, see, e.g., [20]. 

The Cayley transformation maps the imaginary axis into the unit circle. The 
Cayley-exponential function and other new exponential functions presented in this 
paper also have this property. Therefore related T-analogues of trigonometric func- 
tions satisfy exactly the Pythagorean identity: cos 2 t + sin 2 t = 1. 

2 Preliminaries and notation 

In this section we briefly present some preliminaries on time scale calculus, with a 
special stress on standard approaches to exponential, hyperbolic and trigonometric 
functions. 

2.1 Basic notation and definitions 

• Forward jump operator a: a(t) = inf{s € T : s > t) = t a 

• Backward jump operator p: p(t) = sup{s 6 T : s < (} 

• Right-dense points: <r(t) = t. Left-dense points: p(t) = t. 

• Right-scattered points: a(t) > t. Left-scattered points: p(t) < t. 

• Grainincss: p(t) = a(t) — t. 

• Rd- continuous function is, by definition, continuous at right-dense points and 
has a finite limit at left-dense points. 

• Graininess p, is not continuous at points which are left-dense and right-scattered 
but is always rd-continuous. 

• There are two main T-analogues of the i-derivative: 

!{o{t))-f{s) 



Delta derivative f A (t) nm 

s — > t 
s ± a(t) 

Nabla derivative f v (t) := lim 

s t 

S ? p{t) 



a(t) - s 

f(p(t))-m 
Pit) - s 



2.2 Exponential functions 

Delta exponential junction (see [17]), denoted by e a (t,to), is the unique solution of 
the initial value problem 



.r 



A 



a(t)x , x(t ) = 1 , (1) 



where a : T — > C is a given function. Nabla exponential function (see [3]), denoted 
by e a (t, t<j), satisfies: 



y(t) x , x(to) = 1 
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In the continuous case both exponential functions reduces to 

T = M ===> e a (t, t a ) — e a (t, to) — exp / a(s) ds . (2) 

J t 

In particular, 

T = R,a(t) = z e a (t) = e a (t) = e zt . 

In the discrete constant case (T = hZ, a(t) = z G C) we have 

e,(i)=fl + ~) , ^ = ( 1_ f) ' t = nh. 

Another exponential function is related to the so called diamond-alpha derivative 
[25]. However, the associated differential equation is of second order. 

2.3 Hyperbolic and trigonometric functions on time scales 

There exist two different approaches to hyperbolic and trigonometric functions on 
time scales. This ambiguity can be explained as follows. In the continuous case we 
have 

T = R cos x = - — — , sin x 



2 ' 2i 
Unfortunately, the property e~ lx = {e ix )~ l , crucial for some important properties of 
trigonometric functions, does not extend on delta and nabla exponential functions: 

e- a (t,to) / e-\t,t Q ) , e_ a (t,t ) ± e~\t,t Q ) . 

Therefore, trying to extend definitions of trigonometric functions on time scales, one 
has to choose between two natural possibilities: to replace e~ lx either by (e^) -1 or 
by e_; x . Hilger chose the first option [19], i.e., 

ccMt) = e " W+ 2 e " 1(t) , sinh Q (0 = ^W-^W . 

The main advantage of this definition is the identity 

cosh^i) -sinh^(i) = 1 . 

Unfortunately, derivatives cosh^(i), sinh a (t), are not proportional to sinh a (i) 
and cosh Q (i), respectively. Instead, we have 

cosh Q (t) = -f- cosh Q (t) + — — + smh a (t) , 

1 + \xa 1 + \xa 

. , a , , \\iu 2 , a + \\iu 2 

smh^(i) = -f- smh a (t) + — — * cosh Q (t) , 

1 + iia 1 + fia 

where may depend, in general, on t. Another disadvantage is cosh^(i) ^ R 
and sinh^i) ^ iR (for u e R). Therefore, it is difficult to extend this approach on 
trigonomeric functions. However, after some algebraic considerations Hilger arrived 
at surprisingly simple formulae (see [19]) 

Lj(t) = const =$> cos w (i) := cos(wi) , sin u (t) := sin(cjt) . (3) 
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In section 7 the above definition of trigonometric functions on time scales will be 
interpreted in a broader framework of exact discretizations. 

The second (more popular) approach to hyperbolic and trigonometric functions 
has been proposed by Bohner and Petersson [6, 7] 

- / j. \ e a (t) + e- a (t) 
cosh Q (t) = , cos w (i) = cosh iw (t) , 

(4) 

sinh Q (t) = — — - — , sin^fi) = -isinh lw (t) . 

We see that cos w (t) G ffi and sin w (t) G M. for uj(t) G K. What is more, the 
derivatives of hyperbolic and trigonometric functions are good analogues of the 
continuous case: 

smh^ (t) = acosh a (i) , sin^(£) = u;cos w (£) , 

cosh^(t) = asinh Q (£) , cos^(t) = — wsin w (t) . 

Unfortunately, in place of Pythagorean identities we have qualitatively different 
equalities, namely 

cosh^(i) - sinh^(<) = e„ MQ 2(t), cos* (<) + sin* (t) = e^{t) . 

Therefore, functions cos u (£) and sin w (i), defined by (4), are not bounded. 

3 The Cayley-exponential function on T 

We are going to define the Cayley-exponential function as a solution of an ap- 
propriate initial value problem. The definition presented in [10] is equivalent, sec 
Theorem 3.3. 

3.1 Definition of the C-exponential function 

Definition 3.1 The Cayley-exponential (C-exponential) function E a (t, to) satisfies 
the following initial value problem: 

x A (t) = a(t)(x(t)} , x(t ) = l, (5) 

where a is regressive (i.e., [ia ^ ±2) and rd-continuous on T, and 

. , ., x(t)+x(t a ) , , . . x + x* 

(*(*)) = 2 , shOT % : (*> = ■ 

Moreover, we denote E a (t) := E a (t,0). 

In the continuous case the C-exponential function reduces to the usual exponential, 
see (2). In the discrete case we have 

T = hZ, a — const, => E a (t) = 2p — , t = nh . 

V 1 - 27i ta > 

Similar formulas in the discrete case were known a long time ago to Ferrand [16] and 
Duffin [15], and then appeared several times in different contexts [5, 14, 20, 26, 29]. 
They are clearly related to the Cayley transformation, see [20]. 
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3.2 Properties of the C-exponential function 

Theorem 3.2 There is a bijection between Cayley exponential functions and delta 
exponential functions. Namely, E a (t, to) = ep(t, t ), if 

1 + ±/i(i)/3(i) 1 - ±n{t)a(t) 

where pLot ^ ±2 and ///3 ^ — 1. 

Proof: Suppose that x = x(t) satisfies (5). Using x CT — x + px A we obtain 

(1 — -Qfl) X' A = QX . 

Hence, x(t) = ep(t,to), where f3 is given by (6). The assumption pa / ±2 (see Defini- 
tion 3.1) quarantees that j3 exists and + 1 7^ (p-regressivity, see [10]). 

Likewise, suppose that x = xit) satisfies x A = f3x, x(to) = 1. Substituting x = 
x CT — /ix A we get (1 + pP)x A = /?x CT . Adding both equations we obtain 

{2 + p(3)x A = 2f3(x) . 

Therefore, x(t) = e a (t,to), where a is given by (6). The assumption p/3 7^ —1 quarantees 
that aa / -2. □ 

Theorem 3.3 The Cayley- exponential function E a can be expressed as follows 

(7) 



E a (t,t ) := cxp ^ C M ( s )(a(s))As^ 



where Cu,( z ) : = — log f— — f/or A 1 7^ 0), i.e., 2 = — tanh — , arte? Co( z ) 

\i 1 - M 2 



Proof: Using Theorem 3.2 we can apply Hilger's results (the cylinder transformation) 
[17, 19] and the formula (7) follows in a straightforward way. □ 

The following properties of the Cayley-cxponcntial function have been derived 
in [10]. The first formula is clearly related to the Cayley transformation. 

E a (t a ,to) = 7-7 E a (t,to) , E a (t, t ) E a (to,h) = E a (t,ti) , 

1 - ±n{t)a(t) 

E a (t,to) = E a (t,to) , (E a (t,to)) 1 = E- a (t,to) , 

a + (3 



E a (t,t ) Ep(t,t ) = E am p(t,t ) , where a © 



i + ±mW 



Interestingly enough, the last formula is identical with the Lorentz velocity trans- 
formation of special relativity provided that we interpret jj as the speed of light. 
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3.3 Numerical advantages of the C-exponential function 

In the continuous case e a ,e a and E a become identical, see (2). We are going to 
compare their "accuracy" at a right-scattered point t (i.e., t a — t = \x ^ 0), assuming 
a(t a ) = a(t) = a. 

, „ n 1 + 1 , „ 1 , „ 

E a (f, t) = p- = l + an+ -(a/i) 2 + -(a/i) 3 + . . . , 

1 — ±[ia 2 4 

e a (t a , t) = 1 + a/i , e„(f , t) = - — - — = 1 + a^i + (a^) 2 + • ■ ■ 

1 — a/x 

where in the second formula a is evaluated at t a . It means that our assumpation 
(a a = a) is essential. Comparing the above expansions with the continuous case 

exp(a/i) = 1 + afi + \{afif + \ {a^f + ... , 
2 b 

we conclude that E a (t a ,t) is a second-order approximation of cxp(a^), while 
e a {t a and e a (t a ,t) are approximations of the first order. 

4 Cayley-hyperbolic and Cayley-trigonometric func- 
tions on T 

A direct consequence of Definition 3.1 is better theory of hyperbolic and trigono- 
metric fuctions. Caylcy-hyperbolic and Cayley-trigonometric functions are defined 
in a natural way: 

Cosh Q (t) := S«(*)+S- a (t) SmMt) ;= E a (t)-E^(t) 

(8) 

Cos. W :=^ (t)+ /- (t) , SinUt):= E - {t) - E - {t) . 

2 2i 

We have no other choice because the Cayley-exponential function enjoys good prop- 
erties like (E a {t))- 1 = E_ a (t), E & (t) = E a (t) and \E iu (t)\ = 1 for u(t) £ K. 

C-hyperbolic and C-trigonometric functions combine advantages of both Hilger's 
and Bohner-Peterson's approach. 

Theorem 4.1 We assume that a, uj are rd- continuous and afi ^ ±2, ujfi ^ ±2i. 
Then 

Cosh 2 (i) - Sinh 2 (t) = 1 , Cos 2 (t) + Sin 2 (i) = 1 , 

Coshf (t) = a(f) (Sinha(*)) , Sinh^(t) = a(t) (Cosh Q (t)) , 
Cos£(£) = -uj{t) (Sin u (i)) , Sin*(i) - w(t) (Cos u (t)) . 

Proof: Straightforward calculation (see also [10]). □ 
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5 T-analogues of ODE motivated by numerical 
schemes 



Constructing T-analogues of ordinary differential equations (ODEs) one usually 
replaces t-derivatives by delta derivatives [1, 7] or, less frequently, by nabla deriva- 
tives [3]. Thus x A = ax and x v = ax are standard T-analogues of x = ax. The 
results of section 3 suggest, however, that equation x A = a(x) is another (perhaps 
even better) T-analogue of this equation. Certainly T-analogues are not unique. 

In this section we are going to show a correspondence between numerical schemes 
and T-analogues of a given ODE. We consider a general ODE: 

x = f{x, t) , t g T , x(t) e C N , f(x(t), t) e C N , 

and present its T-analogues corresponding to several numerical methods. 



5.1 Euler schemes 

The most popular T-analogues are associated with Euler methods. Delta dynamic 
equations correspond to the forward (explicit) Euler scheme, while nabla dynamic 
equations corcspond to backward (implicit) Euler scheme 

Forward Euler scheme xA (t) = f( x {t)it) > 

Backward Euler scheme x v it) = f(x(t),i) . 

5.2 Trapezoidal rule 

In the autonomous case we have 

Trapezoidal rule x A — | (fix) + fix a j) . 

In the non-autonomous case we can consider at least two different possibilities: 

x A = \ifix,t) + fix"X)) , 

* A = 3 (f{x, t) + f(x-, t) + fix, t-) + f(x°,t°)) 

The first one is related to the standard trapezoidal scheme, the second one is more 
symmetric and has some advantages. Taking f(x,t) = a(t)x we have, respectively, 

x A = (ax) => x a = x , x A = (a)(x) => x a = f-f-f x . 

1 — -^[ia a 1 — jMW 

These equations define next new exponential functions on time scales. The second 
of these functions has better properties because it maps the imaginary axis into the 
unit circle for any function a = ait). 

5.3 Implicit midpoint rule 

In the autonomous case we have 



Implicit midpoint rule x A = / 



x + x° 
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It is not clear how to extend this formula on the non-autonomous case, because in 
the general case \ {t + t a ) T. One may consider, for instance, the following scheme 



5.4 Discrete gradient method 

This is an energy preserving numerical scheme for Hamiltonian systems [22, 23]. 
Here we consider one-dimensional separable systems H(p,q) = T(p) + V(q). The 
equations of motion read 

. _ dT . _ dV 

^ dp 1 ^ dq 

The discrete gradient method yields the following T-analogue of these equations of 
motion: 

a AT A Ay 

= ~~Ap ' = ~Aq (9) 

where the "discrete gradient" is defined as 



Ap y ' P-vp p a -P Aq V ' Q^q q a -Q 

The dynamic system (9) preserves exactly the energy, i.e., H a = H (more detailed 
discussion will be presented elsewhere). We point out that equations (9) differ from 
equations of motion considered in [2], 



5.5 Classical harmonic oscillator equation 

Implicit midpoint, trapezoidal and discrete gradient schemes yield the same T- 
analoguc of the harmonic oscillator equation (g + oj^q = 0): 

AA , 2 // \\ n ,/ \\ I"" + + q 

Q +Uo((q)} = 0, (W» := ^ ■ 

Of course, all these schemes yield different results for other, nonlinear, equations. 

If uj(t) = loq = const, then Cayley-sine and Cay ley-cosine functions satisfy the 
above dynamic equation. Therefore, this T-analogue of the harmonic oscillator has 
bounded solutions. 



6 Pade-exp orient ial functions on T 

The Pade approximation consists in approximation by rational functions of pre- 
scribed order (see, e.g., [4]). Pade approximants Rj,k ar e rational functions 



Rj,k( x ) 



Qk{x) 

which agree with e x (at x = 0) to the highest possible order. Thus, e.g., Rj } k(0) = 1. 
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Definition 6.1 The Pade- exponential function E^ k {t,to) satisfies the dynamic sys- 
tem defined at right dense points by x A = ax and at right scattered points by 

x a = Rjk(afi) x . 

Known exponential functions can be considered as particular Pade exponentials: 
E" (t,to) = e a (f, fo) delta exponential function , x a = (1 + ctfijx , 



Eoi(t,ta) = e a (t,to) nabla exponential function , 

E^^t, to) = E a (t, to) Cay ley- exponential function 
Another special case, x = E^tMolj satisfies equations: 



i 



1 — afx ' 

1 + |ct/j 
1 - \a\i 



x a = 1 + \a^+ jzjati) 2 x __ Q = « x 

l-|a/i+^(a/j) 2 l + ^{ajj) 2 1 - TjCefi + -^(cefi) 2 

Symmetric Pade-exponentials (E k k (t a , to)) map the imaginary axis into the unit 
circle and, therefore, they generate trigonometric functions with "good" properties. 

7 Exact analogues of elementary/special functions 
on T 

Given / : R — > C, we define its exact analogue / : T — > C as / :=/|t, i-e., 

/(f) := /(f) (for t £ T) . 

The path / — > / is obvious and unique, but the inverse way (to find / corresponding 
to a given /) is, in general, neither obvious nor unique. However, in some particular 
cases a natural correspondence / — > / may exist. The constant function / is a 
typical example [10]. The corresponding / is, obviously, constant as well. 

Definition 7.1 We assume a = const 6 C. The exact exponential function on T 

is defined as E^ x (t,t ) := e^*"* '. 

Similarly we can define exact trigonometric and exact hyperbolic functions. Ex- 
act trigonometric functions on T coincide with trigonometric functions on T intro- 
duced by Hilger [19], see also section 2.3. 

Theorem 7.2 The exact exponential function E^ x (t,to) satisifi.es 

x A {t) =a Va(i) (x(t)) , x{t Q ) = l, 
where ip a (t) = 1 at right- dense points and 



, ,s 2 , ault) 
ctfi(t) 2 



at right-scattered points. 



Exact discretizartions of ODEs with constant coefficients were first considerd by 
Potts [28], see also [9, 12, 24]. 
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8 A modification of the g-calculus 



The C-exponential function suggests analogical modifications in the g-calculus. A 
survey of classical results can be found in the textbook [21]. It is curious that 
among several definitions of g-trigonometric functions, used throughout all the pre- 
vious century, none satisfied the Pythagorean identity. Below we present our recent 
results, see [11], enjoying this property. 

New g-exponential function can be defined in two equivalent ways, either as an 
infinite product or an infinite series: 



ex _ TT l + g fc (l-g)f fX V- x_ 

q 11 l-g fe (l-o)| ' q ^ {n}\ ' 

k=0 H V H > 2 n=0 L ' 



where 



{n}\ = {1}{2} . . . {n} , {k} := 1 + g + - ' • + ? 

)[l + q" l ) 



fe-i 



2 1 

Our new g-exponential function can be directly expressed by standard g-exponential 
functions: := eq Eq (for definitions and more details on and see [11, 21]). 

New g-trigonometric functions, motivated by the Cayley transformation, are 
defined by 

Cix c—ix cix i c — ix 

C ■ _ 1 Q n _ 1 1 
OlriqX — , LOSqX — , 

and have the following properties: 

Cos 2 q x + Sin^x = 1 , DqSinqX = (Cos q x) , D q Cos q x = —(Sin q x) , 

where g-derivative is defined by DJ(x) = f<yQ ^ - ^ and (f(x)) = ^ + f<yQ ^ . 

qx — x 2 



9 Dynamic systems on Lie groups 

Two natural generalizations of the Cayley transformation are well known: 

• Lie algebra q — > ( "quadratic" ) Lie group G, 

• anti-Hcrmitcan operators — > unitary operators. 

Here we focus on matrix Lie groups. It is well known that for quadratic Lie groups 
(including all orthogonal, unitary and symplectic groups) we have 

A eg => (I- A)-\I + A) e G . 

Therefore, a natural T-analogue of — <f> = A<f> (here A E g, $ e G) is the dynamic 

at 

equation <i> A = A ($) . The corresponding evolution at right-scattered points is 
expressed by the Cayley transformation: 

„ / + \uA 

A different approach to dynamic systems on Q=SU(2) can be found in [8]. 
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10 Conclusions and future directions 



The most important message of this paper consists in showing that one can define 
many non-equivalent exponential functions on time scales and many non-equivalent 
T-analogues of ordinary differenatial equations. Differential equations have no 
unique 'natural' time scales analogues. It is worthwhile to consider different nu- 
merical schemes in this context, compare [27]. 

Our results suggest numerous developments. We name only few of them: dy- 
namic systems preserving integrals of motion and Lyapunov functions, new de- 
velopments in the g-calculus (e.g., modifications of g-gamma function and of the 
Jackson integral), modifications of g-Laplace and q- Fourier transformations, and 
locally exact T-analogues of elementary functions, compare [13]. 
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